Asymmetric spatial structure of zero modes for birefringent Dirac fermions 
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We study the zero energy modes that arise in an unusual vortex configuration involving both the kinetic 
energy and an appropriate mass term in a model which exhibits birefringent Dirac fermions as its low energy 
excitations. These zero modes only for an appropriate choice of relative vorticities of the mass and kinetic 
energy topological defects. We find the surprising feature that the ratio of the length scales associated with 
states centered on vortex and anti-vortex topological defects can be arbitrarily varied but that fractionalization 
of quantum numbers such as charge is unaffected. We discuss this situation from a symmetry point of view and 
present numerical results for a specific lattice model realization of this scenario. 
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I. INTRODUCTION 

It is well known that the Dirac Hamiltonian with topologi- 
cally non-trivial mass terms allows for zero energy modesJi^ 
If there is a single zero energy mode, then this leads to the 
non-trivial phenomenon of quantum number fractionalization, 
e.g. charge as has been experimentally confirmed in one di- 
mension in polyacetylencj^iiS More recently, there has been 
much interest in zero modes that can arise from topological 
defects in systems whose low energy excitations can be de- 
scribed using Dirac fermionsjil^ Such modes are being stud- 
ied intensely due to the possibiUty of their application in quan- 
tum computationisliiii 

In all of these examples, the zero modes arise as a conse- 
quence of a topological defect in the mass term of the Dirac 
Hamiltonian. We consider here the recently introduced model 
of birefringent fermions^ and introduce a momentum space 
vortex in addition to a standard mass vortex. The low en- 
ergy theory of birefringent fermions consists of four com- 
ponent massless fermions with two separate Fermi velocities 
uo(l ± /?) controlled by the parameter < (3 < 1. Writing the 
low energy theory in Dirac form, the parameter /3 multiplies 
terms in the kinetic energy not present in the regular Dirac 
Hamiltonian. 

Our main result is that when there is an appropriate vortex 
(anti-vortex) in /3 in addition to a mass vortex (anti-vortex) of 
the type introduced in Ref. |^ then the ratio of the characteris- 
tic lengthscales of the zero mode solutions in the presence of 
either a vortex or an anti-vortex may be tuned arbitrarily by (3. 
Specifically, we find that finite /3 breaks the symmetry present 
when /3 ~ between the zero mode solutions in the pres- 
ence of a mass vortex and a mass anti-vortex. The zero mode 
solution in the presence of a vortex becomes more extended 
while with an underlying anti-vortex it becomes more local- 
ized. In the limit (3^0 the model we consider displays the 
same physics as that discussed in Ref. The type of vortex 
we consider here may also be of interest in a variety of sys- 
tems whose low energy excitations can be described as Weyl 
fermions with multiple Fermi velocities, which have recently 
been the focus of a number of publicationsi^— 

This paper is structured as follows. In Sec.|II]we recall the 



model of birefringent Dirac fermions and describe the vortex 
that we are considering. In Sec. |III]we find the zero energy 
modes associated with this vortex, and study a tight-binding 
lattice model to illustrate this physics numerically in Sec. |IV] 
In Sec.[V]we discuss our results and conclude. 



II. MODEL AND VORTEX 

The model of birefringent fermions introduced in Ref. 
has the feature that there are massless fermions near the Dirac 
points, but with two distinct Fermi velocities. This model may 
be obtained as the low-energy theory associated with spinless 
fermions at half filling in a particular tight binding model on 
a square lattice with a four site unit cell illustrated in Fig.[T] 
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FIG. 1. Unit cell of tight binding model with birefringent Dirac 
fermions as low energy excitations.— Hopping parameters are indi- 
cated, with J± = Jo{l ± P). 

The dispersion relation reads as 



±J± 



where J± = Jo(l±/3), with </?<!. This dispersion leads 
to four equivalent Dirac points at the corners of the Brillouin 
zone: K± -t — (±f , ±f )■ Labelling the four sites in the unit 
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cell as A, B, C, and D we can write the low energy theory in 
the form; 

H = Y,i^l[Ek-Hk]^k, (1) 

k 

where V'^ = icAk,CBk,cck,CDk), with cjk a fermionic an- 
nihilation operator for a fermion with momentum k which re- 
sides on sites I ~ A,B, C, or D, and (setting 2 Jg = 1) 

Hk = [i (7071 + «/373) kx + i (7072 + ^,375) ky] . (2) 

We use a non-standard representation of the gamma matri- 
ces in which 70 = 0-3 ® (T3, 71 = 0-2 ® h, 72 = ® (T2, 
73 = -cTi ® ^2, and 75 = 70717273 = -cs (g) cri. The 
matrices 70, 71, 72, 73 and 75 satisfy the Clifford algebra 
7^1 7i^ + li^lfj. = 2(5^,^1^ The representation of the gamma 
matrices is four dimensional, which is the minimal dimen- 
sion for a time-reversal invariant system of spinless Dirac 
fermions in two dimensions on a lattice;^ By way of compar- 
ison, in graphene, the minimal representation is constructed 
with two sublattice degrees of freedom and two inequivalent 
Dirac points. In the present problem, on the other hand, there 
are four equivalent Dirac points, but the unit cell is comprised 
of four lattice points. The two Dirac cones with Fermi ve- 
locities 1 ± /3 suggest the problem may be written as a direct 
sum of two copies of two component massless Dirac fermions 
or Weyl fermions. [This possibility was recently considered 
and it was shown that one cannot have a Weyl fermion in this 
present scenario. -] The realization of Weyl fermions in a two 
dimensional lattice without broken time reversal symmetry is 
also prohibited by Nielsen-Ninomiya theorem.^^^i^ However 
in the limit /3 = 1 there are two flat bands at zero energy 
and a two component Weyl fermion. Such a possibility was 
previously discussed in other contextSii2i2i, _ nevertheless the 
Nielsen-Ninomiya theorem is respected there as well. 

A. Topological Defect 

When /3 = in the birefringent fermion model, there is a 
chiral SU(2) symmetry generated by 73,75, and 735, where 
735 = *7375i but this symmetry is broken when /3 7^ 0,sl If 
we consider the Hamiltonian [Eq. (|2]i] with non-zero (3 and 
introduce a vortex such that /3 ^ in the centre, then chi- 
ral symmetry is restored in the centre of the vortex. Taking 
/3 to be solely a function of the radial co-ordinate r, that van- 
ishes at the origin and has a constant limit as r cxo, then an 
appropriate Hamiltonian for an anti-vortex in (3 reads as 

~P{r) {cos{vpe)^z - sin(t>^6l)75) {-id^) 
-I3{r) (cos(z;^0)75 + sin(t-^0)73) {^idy) (3) 

The vorticity vp is a non-negative integer. One may arrive at 
Eq. (O from Eq. (|2|i by the substitution /3 — > /3(r) and a chiral 
rotation Uc = e*"/^ 07,35 jjj gq Q Without a mass term leading 
to a gap in the spectrum, this vortex in /? does not support a 
normalizable zero mode. 



In order to obtain zero modes, we need a mass term in the 
Dirac-like Hamiltonian to open a gap. An appropriate term 
to consider is the mass term considered by Hou et al^ in the 
context of graphene. In that context, the term corresponds to 
a staggered hopping in a Kekule pattern. We show that this 
term leads to gap for birefringent fermions in AppendixlAl In 
the square lattice problem of Ref. the equivalent term is 
also a staggered hopping, which leads to contribution in the 
low energy Hamiltonian of 

i/„j = -m(r) (cos(w„i6')z7o73 - sin(t;„i6')i7o75) , (4) 

corresponding to a mass anti-vortex with spatial profile given 
by m{r) and vorticity Vm, which is a non-negative integer 
Note that 17073 corresponds to staggered hopping in the x di- 
rection and ^7075 corresponds to staggered hopping in the y 
direction. We will hence work with the Hamiltonian H^.p = 
Hp + H,n- In the limit p ~ Q, the Hamiltonian H^.p takes 
the form of the massive Dirac Hamiltonian originally studied 
by Hou et al^ An index theorem therefore guarantees the exis- 
tence of Vm normalizable zero energy statesM and concomi- 
tant fractionalization of charge when — 1 and /3 = 0^ 
Here we seek to find the zero energy modes in the spectrum 
of H„i^p with /3 ^ Q and study how the anti- vortex in /3 affects 
the zero modes that arise from H„i . 

For calculational convenience, we make use of a unitary 
transformation U ~ U3U2U1, where Ui = I2 (B a'2, U2 = 
{I4 - j(Ti ® I2) , and 1/3 = -i= {I4 + ia2 ® 0-3) to the 

"graphene representation" of the gamma matrices: = 
h ® cTs, 7p = cr3 (g) (J2, ^2 = cTi, 1§ = cTi® 0-2, and 
7^ = (72 cr2 1^ Then the transformed pieces of the Hamil- 
tonian take the form 

Hp = i^Q^i{-idx) - ijo72i-idy) 

-/3{r) cos{vp9)'y3 + sm{vp6)'y5) {~idx) 
-P{r) {cos{vp9)'-f5 + sin(t;^6')73) {-idy), (5) 

and 

Hm = -m{r) (- cos(i;„i6')i7o73 - sin(i;„i6')i7o75) • (6) 

We assume that the spatial profiles of the defects in /3 and m 
are that as r — > 0, /3{r) — > and m(r) 0, and in the 
large r limit (far from the core of the vortex), /3(r) — > /3o and 
m{r) — >■ niQ, where /3o and niQ are constants. It might be pos- 
sible to realize such a vortex configuration experimentally by 
implementing a pattern of hopping integrals such as described 
in Sec. HV] for cold atoms in an optical lattice^**^' or through 
other synthetically constructed systems with Dirac fermionic 
excitations 4^ 



III. ZERO ENERGY MODES 

We now search for zero energy modes that satisfy 
Hm,p'^ = 0. The existence of a unitary operator, 70, such 
that {H,n.p:lo] ~ 0, ensures the spectral symmetry of the 
energy eigenstates of Hm,p- Moreover, one may find an anti- 
unitary operator M = UK, where U is unitary and K is 
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the complex conjugation operator, which anticommutes with 
the Hamiltonianjii Noting that {17071, 75, *7o73} real and 
{«7o72,73,«7o75} are imaginary one finds U = -^7273 = 
(Ji ® (73, in the "graphene representation". A zero energy 
mode of Hp^„i must also, therefore be an eigenstate of A/, 
which for the state = (-01, ip2,ip3,ipi) leads to the con- 
straint -01 = ±V'3' ^2 = TV'I- The symmetry of the energy 
spectrum about zero implies that the zero energy mode should 
be robust against any weak local perturbation. 

For further ease of calculation we redefine the components 
of as eT7ov]/ so that the eigenvalue equations for the 
zero energy mode take the form 



dzi^2 - (/3(r)e-*"^^aj - m(r)e- 
9201 - (/3(r)e-™^'*aj + m(r)e- 
9,-04 - (/3(r)e™^«a, - m(r) 
9,^3- (/3We™^'a,+m(r) 



04 = 0, (7) 

03 = 0, (8) 

^2 = 0, (9) 

^1 = 0, (10) 



where = d^-idy = e"'^ [dr - ^de) and = d^+idy = 
e*^ [dr + ^dg). ^ must be an eigenvector of M as well as H, 
and choosing to have eigenvalue -1 implies 03 = — 0j[' and 
04 = 1/^2 (we will discuss the choice that M has eigenvalue 
+1 below). Under this constraint Eqs. (|7]i - (fTol i lead to only 
two independent equations: 

- {I3{r)e-'''^'d-, - m(r)e-"'"^) ^2 = 0, (H) 
d-zi^i + (/3(r)e-"'''^9. + m(r)e-™'"^) 0-* = 0. (12) 

To solve these equations we make use of an ansatz introduced 
in a different context by Ghaemi and Wilczek,"*^^ focusing first 
on Eq. (fTTI) . We make the ansatz 



02(r,0) = e^'V2(r) + e™>4(r), 



(13) 



where 02 and 04 are real. In order to have a consistent so- 
lution, we must have that Vm = — 1. There is still some 
freedom in the choice of I and n, depending on which terms 
are grouped together from Eq. ( fTTI ) after the use of the ansatz 
in Eq. ( fT3] l. The choice that allows either I = n 01 1 ^ n and 
guarantees that the solution is single- valued is: 



dr 
dr 



<f>2 - m 

04 - /3(r) 



dr 
dr 



04 + m(r)04 0, (14) 
02+m(r)02 =0, (15) 



in which case we have the condition I + n = 2 — vp. We 
now consider the asymptotic behavior of the solution of this 
equation at large and small r. In the large r limit j3{r) — > /3o 
and 771(7-) — > ?77o, and we can ignore 1/r terms. This leads to 
the solution as r ^> 00 



where 



^2(0 =04(r) = Ae-^-^ 



rriQ 
l~/?o' 



(16) 



is the inverse of the characteristic length scale for the zero 
mode in the presence of an anti-vortex and ^ is a normaliza- 
tion constant. At small r, I5{r) — > and 777(r) — > 0, in which 
case 



02 (r) = a2 



04 (r) = 04 r 



(17) 



In order that solution be normalizable near the origin, we re- 
quire ^, 77 < 0, which implies U/3 > 2. 

We now consider the two cases vp odd and even sepa- 
rately. If vp is even, i.e. vp = 2p for some integer p, 
then to satisfy Eqs. (fT4l i and ( fTSI ) subject to I, n < and 
I + n = 2 — vp = 2(1 — p) requires p > 1. There will be 
a solution with I = n and (2p — 2) solutions with I ^ n. As 
far as the original equation we were trying to solve, Eq. (fTTI) . 
is concerned, solutions with I ^ n are equivalent, so there 
are in fact only p — \ solutions of Eq. ( fTTI ) with I ^ n, and a 
total of p zero mode solutions when vp = 2p (which implies 
v,n = 2p — 1). When vp is odd, we can apply a similar anal- 
ysis and find that for vp — 2p + 1 (i.e. u,n = 2p), then there 
are p zero mode solutions. Hence to have a single zero mode 
with M eigenvalue —1, we can either have vp — 2, Vm — 1, 
or Vp = 3, Vm = 2. 

In our solution above we assumed that the eigenvalue of Af 
is —1. If, on the other hand, we try to obtain a solution which 
is an eigenvector of M with eigenvalue +1, then we find that 
we must have vp > 2 and Vm = vp — 1, as before, and that 
there is no normalizable solution when vp = 2. However, 
there can be normalizable solutions when vp > 2, and in par- 
ticular there is one normalizable solution when vp = 3, so 
that there are a total of 2 normalizable zero modes (as would 
be expected from Vm = 2 by the usual index theorem^). We 
consider the situation with vp ~ 2, Vm = ^ numerically in 
Sec. |IV] It should be noted that the fact that the solution we 
obtained from Eqs. ( fTTI ) and ( fTSl i has eigenvalue —1 for M 
is not significant. We could have equally well found a nor- 
malizable solution with eigenvalue +1 and no normalizable 
solution with —1 when vp = 2 had we made a different redefi- 
nition of the fields before Eqs. (|7]i to ([TOll, e.g. * e^^'^o*. 

We can also see that the mass anti-vortex is required in or- 
der to have normalizable solutions. If we set 7770 = 0, then 
we have an anti-vortex defect in the kinetic energy alone. At 
large distance, the solutions asymptote to a constant value as 



02 (r) = 04(0 ^ A, 



(18) 



leading to a non-normalizable solution. Note that additionally, 
when 7770 ~ there is no gap in the spectrum. When vp > 0, 
there are no non-zero normalizable solutions for 

Next we consider vortex configurations in /3 and m, corre- 
sponding to negative integer values of vp and Vm, which will 
illustrate how the /3 term modifies the usual solution to the 
mass vortex;^ In this case, 0^2 has no normalizable solutions. 
Nevertheless, one can obtain the zero energy mode by using 
the ansatz 



0i(r,(?) =e*"'0i(r)+e-''03(r), 



(19) 
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and ■03 = —ipl, associated with the —1 eigenvalue of M. The 
coupled differential equations for the zero mode then read as 



03 + m{r)h = 0, (20) 
01 +m(r)0i = 0. (21) 



Even though in the vicinity of the origin 0i and 03 behave 
identically to 02 and 04 respectively in Eq. (fTTb . at large dis- 
tances (r — !• oo) 
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(22) 



where 



mo 

+ l + /?o 

is the inverse of the characteristic lengthscale for the zero 
mode in the presence of a vortex, and i? is a normalization 
constant. The lengthscale for the decay of the zero mode dif- 
fers in this case from the lengthscale we found for the anti- 
vortex solution. The norm of the zero modes with +1 eigen- 
value of M in this situation grow with the system size when 

V/^ = -2, Vm = -1- 

In the limit (3q — 0, the characteristic length scales for vor- 
tex and anti-vortex are the same, since the Hamiltonian asso- 
ciated with those two distinct topological defects are unitarily 
equivalent to each other. Namely, the vortex Hamiltonian may 
be obtained from the anti-vortex one via a unitary rotation by 
75 in Eq. (|6]l. However, this unitary equivalence breaks down 
for a birefringent Dirac Hamiltonian, leading to the distinct 
length scales found in Eqs. ( fTSI ) and i2% . 

It is useful to present the zero modes in terms of the orig- 
inal spinor components ipk- In the presence of an anti- vortex 
the zero mode reads as ipQ = (0, -02, 0, '02)- After unitary ro- 
tation hyU = UlU^Ul one finds that the amplitudes of the 
zero energy mode with an underlying anti-vortex is finite only 
on the B and C sublattices. When there is an underlying vor- 
tex, the zero modes acquire a finite expectation value on the A 
and D sublattices. This can also be seen by noting that both 
the vortex and anti-vortex Hamiltonians anticommute with 70. 
Therefore the zero energy subspace "Hq is invariant under 70 
and it acts like an identity matrix in Hq. Therefore all the zero 
energy states must be an eigenstate of 70 with eigenvalue +1 
or —1. Since 70 = (73 ® (T3 = Diag(l, —1, —1, 1), the am- 
plitude of the zero modes can be finite either on the B and 
C sublattices (anti-vortex) or A and D sublattices (vortex) in 
accordance with our explicit calculation. 



we have obtained in the continuum limit carry over to the lat- 
tice. To do this we considered the lattice model introduced in 
Ref. I23 and introduced a Z4 vortex (anti-vortex) with vortic- 
ity 2 in /3 and a Z4 vortex (anti-vortex) in to, as illustrated in 
Fig.|2]for antivortices in (3 and to centered on an A site. 
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FIG. 2. Contributions to the hopping integrals from a) Anti-vortex in 
fi term and b) Anti- vortex in m term. The color coding in a) is blue 
for J_, red for J+, green for — J+ and white for — J_. The color 
coding in b) is red for m and blue for — m (note that not all links 
have a contribution to the hopping in case b)). We have suppressed 
the radial dependence of the hopping parameters for simplicity. The 
central site is an A site. 



IV. NUMERICS 

The results we obtained regarding the zero modes in Sec. Hill 
are for a continuum theory. Previous studies^i^ have found 
that zero modes that are predicted from a continuum calcula- 
tion are present in the spectrum of appropriate tight-binding 
models on a lattice. Given the complicated form of the vortex 
configuration we considered here, we check that the results 



We diagonalized the Hamiltonian on an L x L lattice and 
considered topological defects centered on A, B, C and D 
sites. We chose a vortex (anti-vortex) in both f3 and m with 
step-function spatial profile /(r) oc 9{r — ro), with ro = 3 
lattice spacings. Our results were qualitatively similar with 
topological defects centred on different sites, and we mainly 
display data for configurations of the type shown in Fig. |2]in 
which there is a topological defect centered on an A site. The 
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spectrum as a function of system size is shown in Fig. [3] for 
mo = 0.25, /3q ~ 0.25 for a topological defect centered on 
an A site. We see that there are two states which converge 
to zero energy, and a series of non-zero energy states in the 
gap - the continuum levels are visible at the top and bottom of 
the energy range. We can also see that there are small differ- 
ences in the spectrum when there are vortex and anti-vortex 
configurations of the hopping. 




centered on site A. 



Even though there is only one vortex (anti-vortex) on the 
lattice, there are two zero modes in the numerical spectrum 
as is expected for a finite size system with open boundary 
conditions4i The probability density of these states is local- 
ized at the centre of the lattice on the vortex (anti-vortex) and 
at the edge of the system. The integrated charge density av- 
eraged over all four possible locations of the topological de- 
fect (A, B, C or D) for the zero energy modes is illustrated in 
Fig. |4] We compare the cases mo = 0.25 and = and 
jTio = 0.25 and /?o = 0.75, and it is quite evident that there 
is precisely half a charge localized around the centre of the 
system. The profiles also are in qualitative agreement with 
the spatial dependences derived in Sec.|IIIl The antivortex so- 
lution with Po = 0.75 is considerably more localised than 
the vortex solution with the same /3o and the solution when 
/3o = 0, both at small r and at larger r > 3ro (see the inset to 
Fig. HI. 

A feature of our analytical results is that the vortex and an- 
tivortex solutions have support on different sublattices. This 
is also borne out in our numerical results. We show the charge 
density for a vortex with mg = 0.25 and /3o = 0.75 in Fig.|5] 
in which there is support for the state only on A and D sites, 
as deduced in Sec. |III] We also confirmed that the state in the 
presence of an anti-vortex only has support on B and C sites. 




r 



FIG. 4. Integrated charge density near the centre of a 102 x 102 size 
system showing the spatial profile of the vortex and anti-vortex so- 
lutions with mo = 0.25 and /3o = 0, and the spatial profile of the 
vortex and anti-vortex solutions when mo — 0.25 and /3o = 0.75, 
averaged over defects located on A, B, C, and D sites. 




FIG. 5. Charge density near the centre of a 102x 102 size system 
centered on an A site showing the spatial profile of the vortex so- 
lutions when mo = 0.25 and /3o ~ 0.75 (the core of the vortex is 
shaded). Note that there is support for the state only on A and D 
sites. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the zero modes of birefringent 
Dirac fermions arising when there is an unusual vortex con- 
figuration which involves the kinetic energy as well as a mass 
term. This situation differs from the usual one in which there 
is fractionalization of non-interacting fermions in one or two 
dimensions due to a toplogical defect in some mass term,"*^ In 
order to have a normalizable single-valued solution, there is a 
constraint on the allowed vorticity of the vortex (anti-vortex) 
in (3 in that it winds one more time (in the same sense) as the 
vortex (anti-vortex) in m. 

The effect of the topological defect in /3 that we consider 
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is to differentiate the spatial profiles of the zero modes asso- 
ciated with vortex and anti- vortex zero modes. The presence 
of the toplogical defect in the kinetic energy in addition to 
the mass leads to the unusual situation in which the length 
scale associated with the localized state at a vortex differs 
from the length scale associated with a state centered on an 
anti-vortex: the ratio of the two characteristic length scales 
is (1 + /3o)/(l — f3o) and so can be made arbitrarily large. 
The role of the vortex and the anti-vortex can be interchanged 
by changing the sign of /?. As we emphasised in Sec. |III] 
the introduction of (3 breaks the chiral symmetry that relates 
the vortex and antivortex Hamiltonians, allowing the possi- 
bility of differing zero mode solutions. It should be noted 
that the topological defect in (3 does not affect the fractional- 
ization associated with the topological defects, only the rele- 
vant lengthscales. We demonstrated this feature qualitatively 
through diagonalizing a lattice model and comparing the inte- 
grated charge density for the cases of /3o = and finite /3o. 



The results we have obtained here have interest in a broader 
context than the particular problem we studied. Our work 
gives an example of a physical property that systems which 



recently discussed by one of us in the context of graphene 
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,26-30 



can 



exhibit birefringent massless fermionic excitation; 
have that are unavailable for regular Dirac fermions (an other 
such example is birefringent Klein tunnelling^^). The unusual 
occurrence of zero modes with different lengthscales (but not 
with associated fractionalization) was noted by two of us in 
Ref. in the context of zero modes for birefringent Dirac 
fermions in the presence of a domain wall in a 70 mass term, 
and may be a generic feature of zero modes in birefringent 
Dirac systems. A situation in which there are two zero modes 
(hence there is no fractionalization) associated with a topo- 
logical defect, but each has a differing lengthscale was also 
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Appendix A: Dispersion of birefringent Dirac fermions in the 
presence of a 17073 mass term. 

In order for there to be localized zero modes in the presence 
of a topological defect, we need to be certain that the addi- 
tional term in the Hamiltonian leads to a gap in the spectrum 
in the absence of such a defect. To confirm this, we calculate 
the spectrum for the following Hamiltonian: 

(Al) 

After a short calculation one may determine that the eigenval- 
ues are 



efe = ± 



\ 



(l+/32)|fc|2+77i2±2^|/c| 




Note that there is always a gap of 2m at fc = 0, and that 
the minimum energy can occur at a finite value of k. For 
appropriate choices of /? and m there can be a gap, and the 
choices we make for our numerical calculations in Sec. |IV] 
are such that a gap always exists. An analogous calculation 
with an 17075 mass will yield a result such as this with mky 
replaced by mkx- 
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